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WHAT IS CLAIMED IS: 

1- . An internal exhaust gas recirculation amount estimation 
system of an internal combustion engine comprising: 

an exhaust valve closure timing in- cylinder temperature 
5 calculation section that calculates an in-cylinder 
temperature at an exhaust valve closure timing; 

an exhaust valve closure timing in-cylinder pressure 
calculation section that calculates an in-cylinder pressure 
at the exhaust valve closure timing; 
10 a gas constant calculation section that calculates a gas 

constant corresponding to a change in a composition of 
exhaust gas, based on an air-fuel mixture ratio; 

an exhaust valve closure timing in-cylinder residual gas 
amount calculation section that calculates an in-cylinder 
15 residual gas amount at the exhaust valve closure timing, 
based on at least the in-cylinder temperature, the in- 
cylinder pressure, and the gas constant; 

an valve overlap period blow-back gas amount calculation 
section that calculates a valve overlap period blow-back gas 
20 amount defined as a quantity of gas flow from one of intake 
and exhaust ports via a combustion chamber to the other port 
during a valve overlap period during which both of intake 
and exhaust valves are open together, the gas flow being 
created by a pressure difference between the intake and 
25 exhaust ports; and 

an internal EGR amount calculation section that 
calculates an internal exhaust gas recirculation amount 
based on the exhaust valve closure timing in-cylinder 
residual gas amount and the valve overlap period blow-back 
30 gas amount. 

2. The internal exhaust gas recirculation amount 
estimation system as claimed in claim 1, wherein: 
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the internal exhaust gas recirculation amount is 
arithmetically calculated by adding the valve overlap period 
blow-back gas amount to the exhaust valve closure timing in- 
cylinder residual gas amount, from a predetermined physical 
5 expression MRES = MRESCYL + MRESOL, where MRES denotes the 
internal exhaust gas recirculation amount, MRESOL denotes 
the valve overlap period blow-back gas amount, and MRESCYL 
denotes the exhaust valve closure timing in- cylinder 
residual gas amount. 

10 

3. The internal exhaust gas recirculation amount 
estimation system as claimed in claim 2, wherein: 

a value of the valve overlap period blow-back gas amount 
(MRESOL) is defined as a positive value when a gas flows 

15 from the exhaust port toward the combustion chamber and then 
blows from the combustion chamber back to the intake port 
during the valve overlap period, and the value of the valve 
overlap period blow-back gas amount is defined as a negative 
value when a gas flows from the intake port toward the 

20 combustion chamber and then blows from the combustion 

chamber through the exhaust port during the valve overlap 
period. 

4. The internal exhaust gas recirculation amount 
25 estimation system as claimed in claim 1, wherein; 

the exhaust valve closure timing in- cylinder residual gas 
amount calculation section comprises: 

(a) an exhaust valve closure timing in-cylinder 
volumetric capacity calculation section that calculates an 
30 in-cylinder volumetric capacity at the exhaust valve closure 
timing; 
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(b) the exhaust valve closure timing in-cylinder 
temperature calculation section that calculates the in- 
cylinder temperature at the exhaust valve closure timing; 

(c) the exhaust valve closure timing in-cylinder 

5 pressure calculation section that calculates the in-cylinder 
pressure at the exhaust valve closure timing; and 

(d) the gas constant calculation section that 
calculates the gas constant corresponding to the change in 
the composition of exhaust gas, based on the air-fuel 

10 mixture ratio; 

the exhaust valve closure timing in-cylinder residual gas 
amount calculation section calculates the exhaust valve 
closure timing in-cylinder residual gas amount, based on the 
in-cylinder volumetric capacity, the in-cylinder temperature, 

15 the in-cylinder pressure, and the gas constant. 



5. The internal exhaust gas recirculation amount 
estimation system as claimed in claim 4, wherein: 

the exhaust valve closure timing in-cylinder residual gas 
20 amount is arithmetically calculated based on the in-cylinder 
volumetric capacity, the in-cylinder temperature, the in- 
cylinder pressure, and the gas constant, from a 
predetermined physical expression MRESCYL = 

(PEVCxVEVC) /(REXxTEVC) , where PEVC denotes the in-cylinder 
25 pressure, VEVC denotes the in-cylinder volumetric capacity, 
REX denotes the gas constant, and TEVC denotes the in- 
cylinder temperature . 



6. The internal exhaust gas recirculation amount 
30 estimation system as claimed in claim 4, wherein: 

the exhaust valve closure timing in-cylinder volumetric 
capacity calculation section geometrically determines the 
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in-cylinder volumetric capacity based on a position of a 
reciprocating piston at the exhaust valve closure timing. 

7. The internal exhaust gas recirculation amount 
5 estimation system as claimed in claim 4, wherein: 

in the internal combustion engine employing both a 
variable valve timing mechanism that variably adjusts at 
least the exhaust valve closure timing and a variable 
compression ratio mechanism that variably adjusts a 

10 compression ratio by changing at least one of a top dead 
center and a bottom dead center of a piston stroke 
characteristic of a reciprocating piston, the exhaust valve 
closure timing in-cylinder volumetric capacity calculation 
section calculates the in-cylinder volumetric capacity based 

15 on a variation of the exhaust valve closure timing and the 
variation of the compression ratio. 

8 . The internal exhaust gas recirculation amount 
estimation system as claimed in claim 1, wherein: 

20 the gas constant calculation section calculates the gas 

constant corresponding to the change in the composition of 
exhaust gas, based on a target combustion equivalent ratio. 

9 . The internal exhaust gas recirculation amount 
25 estimation system as claimed in claim 1, wherein: 

the valve overlap period blow-back gas amount calculation 
section comprises : 

(a) the exhaust valve closure timing in-cylinder 
temperature calculation section that calculates the in- 

30 cylinder temperature at the exhaust valve closure timing; 

(b) the exhaust valve closure timing in-cylinder 
pressure calculation section that calculates the in-cylinder 
pressure at the exhaust valve closure timing; 
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(c) the gas constant calculation section that 
calculates the gas constant corresponding to the change in 
the composition of exhaust gas , based on the air-fuel 
mixture ratio; 

5 (d) an intake pressure calculation section that 

calculates an intake pressure; 

(e) a gas ratio-of -specif ic-heat calculation section 
that calculates a ratio of specific heat of gas, 
corresponding to the change in the composition of gas 

10 flowing between the intake and exhaust ports via the 
combustion chamber during the valve overlap period; 

(f ) a valve overlap period integrated effective area 
calculation section that calculates a valve overlap period 
integrated effective area for an opening area of the intake 

15 valve and an opening area of the exhaust valve during the 
valve overlap period; 

(g) an engine speed calculation section that calculates 
an engine speed; and 

(h) a supercharging-and-choking phenomena determination 
20 section that determines the presence or absence of a 

supercharging phenomenon regarding the gas flow between the 
intake and exhaust ports via the combustion chamber during 
the valve overlap period, the gas flow being created by the 
pressure difference between the intake and exhaust ports, 

25 and the presence or absence of a choking phenomenon 

regarding the gas flow between the intake and exhaust ports 
via the combustion chamber during the valve overlap period; 

the valve overlap period blow-back gas amount calculation 
section calculates the valve overlap period blow-back gas 

30 amount, based on the in-cylinder temperature, the in- 

cylinder pressure, the gas constant, the intake pressure, 
the ratio of specific heat, the valve overlap period 
integrated effective area, and the engine speed, and a 
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decision result of the supercharging-and-choking phenomena 
determination section . 

10. The internal exhaust gas recirculation amount 
5 estimation system as claimed in claim 9, wherein: 

the valve overlap period blow-back gas amount is 
calculated based on the engine speed, the valve overlap 
period integrated effective area, and a temporary valve 
overlap period blow-back gas amount determined based on the 

10 supercharging-and-choking phenomena decision result, from a 
predetermined expression: 

MRESOL = (MRESOLtmpxASUMOLx60)/(NRPMx360) 
where MRESOLtmp denotes the temporary valve overlap period 
blow-back gas amount, ASUMOL denotes the valve overlap 

15 period integrated effective area, and NRPM denotes the 
engine speed. 

11. The internal exhaust gas recirculation amount 
estimation system as claimed in claim 10, wherein: 

20 the supercharging-and-choking phenomena determination 

section comprises an intake-pressure-to-exhaust-pressure 
ratio calculation section that calculates a pressure ratio 
of intake pressure to exhaust pressure, based on the intake 
pressure and the in-cylinder pressure from a predetermined 

25 expression PINBYEX=PIN/PEVC , where PINBYEX denotes the 

pressure ratio, PIN denotes the intake pressure, and PEVC 
denotes the in-cylinder pressure; 

the supercharging-and-choking phenomena determination 
section determines that the supercharging phenomenon is 

30 present when the pressure ratio is greater than a 

predetermined threshold value, and determines that the 
supercharging phenomenon is absent when the pressure ratio 
is less than or equal to the predetermined threshold value; 
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the valve overlap period blow-back gas amount calculation 
section calculates the valve overlap period blow-back gas 
amount by using the ratio of specific heat of exhaust gas, 
corresponding to the change in the composition of exhaust 
5 gas, when the supercharging phenomenon is absent; and 

the valve overlap period blow-back gas amount calculation 
section calculates the valve overlap period blow-back gas 
amount by using the ratio of specific heat of air-fuel 
mixture, corresponding to the change in the composition of 
10 air-fuel mixture, when the supercharging phenomenon is 
present . 

12. The internal exhaust gas recirculation amount 
estimation system as claimed in claim 9, wherein: 

15 the valve overlap period integrated effective area 

calculation section calculates a valve overlap based on an 
intake valve open timing and the exhaust valve closure 
timing, and calculates the valve overlap period integrated 
effective area based on the valve overlap. 

20 

13. The internal exhaust gas recirculation amount 
estimation system as claimed in claim 9, wherein: 

the valve overlap period integrated effective area 
calculation section sets a smaller one of the opening area 

2 5 of the intake valve and the opening area of exhaust valve as 
an effective valve opening area, and calculates the valve 
overlap period integrated effective area as an integrated 
value of the smaller opening area per unit time from a point 
of time corresponding to a beginning-of -valve -overlap crank 

30 angle to a point of time corresponding to an end-of -valve- 
overlap crank angle . 
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14. The internal exhaust gas recirculation amount 
estimation system as claimed in claim 1, further comprising: 

a variable valve timing mechanism enabling at least one 
of an intake valve open timing and the exhaust valve closure 
5 timing to be varied; and 

a valve timing calculation section that estimates the 
valve timing based on a valve-timing variation adjusted by 
the variable valve timing mechanism and affected by a valve 
temperature change and a motion-transmitting valve-train 
10 component parts wear; 

the valve timing calculation section comprising: 

(a) a valve clearance estimation section that estimates 
a valve clearance; and 

(b) a valve timing compensation section that 
15 compensates for the valve timing based on the valve 

clearance estimated. 

15. The internal exhaust gas recirculation amount 
estimation system as claimed in claim 14, wherein: 

20 the valve clearance estimation section comprises a valve 

temperature calculation section that calculates a valve 
temperature; and 

the valve clearance estimation section estimates the 
valve clearance based on the valve temperature. 

25 

16. The internal exhaust gas recirculation amount 
estimation system as claimed in claim 15, wherein: 

the valve temperature calculation section comprises an 
engine load calculation section that calculates an engine 
30 load; and 

the valve clearance estimation section estimates the 
valve clearance based on the engine load. 
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17. The internal exhaust gas recirculation amount 
estimation system as claimed in claim 14, wherein: 
the valve timing calculation section comprising: 

(a) an exhaust valve clearance estimation section that 
5 estimates an exhaust valve clearance; 

(b) an exhaust valve closure timing compensation 
section that compensates for the exhaust valve closure 
timing based on the exhaust valve clearance estimated; 

(c) an intake valve clearance estimation section that 
10 estimates an intake valve clearance; and 

(d) an intake valve open timing compensation section 
that compensates for the intake valve open timing based on 
the intake valve clearance estimated; 

15 18, The internal exhaust gas recirculation amount 
estimation system as claimed in claim 14, wherein: 
the valve timing compensation section comprises: 

(a) a valve-timing real value calculation section that 
calculates a valve-timing real value based on a frequency of 

20 oscillatory seating motion of each of the intake and exhaust 
valves under a predetermined condition that the engine is 
conditioned in a valve timing detectable engine operating 
state suited to detect the valve- timing real value; 

(b) a valve- timing estimate calculation section that 
25 calculates, under the predetermined condition, a valve- 
timing estimate of the valve timing compensated for based on 
the valve clearance estimated and temperature- compensated 
for based on a change in a valve temperature; 

(c) a valve timing error calculation section that 
30 calculates a valve timing error between the valve- timing 

real value and the valve- timing estimate; and 

(d) a wear dependent valve- timing-correction learned 
value calculation section that calculates a valve-timing- 
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correction learned value based on the valve timing error; 
and 

wherein the valve timing is further compensated for based 
on the valve- timing-correction learned value. 

5 

19. The internal exhaust gas recirculation amount 
estimation system as claimed in claim 18, further 
comprising: 

a knock sensor that detects the frequency of oscillatory 
10 seating motion of each of the intake and exhaust valves. 



20. The internal exhaust gas recirculation amount 
estimation system as claimed in claim 11, wherein: 

the supercharging-and-choking phenomena determination 

15 section comprises a choking-phenomenon decision threshold 

value calculation section that calculates a minimum choking- 
phenomenon decision threshold value SLCHOKEL and a maximum 
choking-phenomenon decision threshold value SLCHOKEH based 
on the ratio of specific heat of gas, from predetermined 

20 expressions: 

SLCHOKEL = {2/(SHEATR+l)} {SHEATR/(SHEATR - 1)} 
SLCHOKEH = {2/(SHEATR+l)} { - SHEATR/(SHEATR - 1)} 
where SHEATR denotes the ratio of specific heat of exhaust 
gas, corresponding to the change in the composition of 

25 exhaust gas, obtained when the supercharging phenomenon is 
absent, and the ratio of specific heat of exhaust gas is 
replaced with the ratio of specific heat of air-fuel mixture, 
obtained when the supercharging phenomenon is present; 
the supercharging-and-choking phenomena determination 

30 section determines that the choking phenomenon is present 
when either of a condition that the pressure ratio is less 
than the minimum choking-phenomenon decision threshold value 
SLCHOKEL and a condition that the pressure ratio is greater 
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than the maximum choking-phenomenon decision threshold value 
SLCHOKEH is satisfied, and determines that the choking 
phenomenon is absent when the pressure ratio is greater than 
or equal to the minimum choking-phenomenon decision 
5 threshold value SLCHOKEL and the pressure ratio is less than 
or equal to the maximum choking-phenomenon decision 
threshold value SLCHOKEH; and 

the valve overlap period blow-back gas amount calculation 
section comprises an unsupercharged and unchoked state, 
10 valve overlap period blow-back gas amount calculation 

section that calculates an unsupercharged and unchoked state, 
temporary valve overlap period blow-back gas amount from a 
predetermined expression: 

MRESOLtmp = 4l xPEVCxMRSOLDxMRSOLP 
15 where PEVC denotes the in-cylinder pressure, MRSOLD denotes 

a density term and MRSOLP denotes a pressure difference term, 
and the density term MRSOLD and the pressure difference term 
MRSOLP are represented by predetermined expressions: 
MRSOLD = SQRT{l/(REXxTEVC)} 
20 MRSOLP = SQRT [ SHEATR/ ( SHEATR- 1 ) 

x { P I NBYEX ( 2 1 SHEATR ) -PINBYEX ( ( SHEATR+ 1 } 1 SHEATR ) } ] 
where REX denotes the gas constant, TEVC denotes the in- 
cylinder temperature, SQRT denotes a coefficient for the gas 
constant REX and the in-cylinder temperature TEVC, SHEATR 
25 denotes the ratio of specific heat, and PINBYEX denotes the 
pressure ratio. 

21. The internal exhaust gas recirculation amount 
estimation system as claimed in claim 20, wherein; 
30 the valve overlap period blow-back gas amount calculation 

section comprises an unsupercharged and choked state, valve 
overlap period blow-back gas amount calculation section that 
calculates an unsupercharged and choked state, temporary 
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valve overlap period blow-back gas amount from a 
predetermined expression: 

MRESOLtmp = PEVCxMRSOLDxMRSOLPC 
where PEVC denotes the in-cylinder pressure, MRSOLD denotes 
the density term and MRSOLPC denotes a pressure difference 
term, and the density term MRSOLD and the pressure 
difference term MRSOLPC are represented by predetermined 
expressions : 

MRSOLD = SQRT{l/(REXxTEVC)} 



where REX denotes the gas constant, TEVC denotes the in- 
cylinder temperature, SQRT denotes the coefficient for the 
gas constant REX and the in-cylinder temperature TEVC, and 
SHEATR denotes the ratio of specific heat. 

22. The internal exhaust gas recirculation amount 
estimation system as claimed in claim 21, wherein: 

the valve overlap period blow-back gas amount calculation 
section comprises a supercharged and unchoked state, valve 
overlap period blow-back gas amount calculation section that 
calculates a supercharged and unchoked state, temporary 
valve overlap period blow-back gas amount from a 
predetermined expression: 



where PIN denotes the intake pressure, RIN denotes a gas 
constant of intake air, TIN denotes a Kelvin absolute 
temperature of the intake air, and MRSOLPT denotes a 
pressure difference term represented by a predetermined 
expression : 

MRSOLPT = SQRT [ SHEATR / ( SHEATR - 1 ) 



MRSOLPC = SQRT [ SHE ATRx { 2 / ( SHEATR+ 1 ) } 



( ( SHEATR+ 1 ) / ( SHEATR - 1 ) } 



MRESOLtmp 




x{PINBYEX ( " 2/SHEATR, -PINBYEX <-(sheatr + i)/sheatr) } ] 
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where SQRT denotes the coefficient, SHEATR denotes the ratio 
of specific heat, and PINBYEX denotes the pressure ratio. 

23. The internal exhaust gas recirculation amount 
5 estimation system as claimed in claim 22, wherein: 

the valve overlap period blow-back gas amount calculation 
section comprises a supercharged and choked state, valve 
overlap period blow-back gas amount calculation section that 
calculates a supercharged and choked state, temporary valve 
10 overlap period blow-back gas amount from a predetermined 
expression : 

MRESOLtmp = - J xPINxMRSOLPC 

\RINxTIN 

where PIN denotes the intake pressure, RIN denotes the gas 
constant of the intake air, TIN denotes the Kelvin absolute 
15 temperature of the intake air, and MRSOLPC denotes a 

pressure difference term represented by the predetermined 
expression : 

MRSOLPC = SQRT[SHEATRx{2/(SHEATR+l) }(<sheatr+i)/<sheatr-i)> j 
where SQRT denotes the coefficient and SHEATR denotes the 
20 ratio of specific heat. 

24. An internal exhaust gas recirculation amount estimation 
system of an internal combustion engine comprising: 

exhaust valve closure timing in-cylinder temperature 
25 calculation means for calculating an in-cylinder temperature 
at an exhaust valve closure timing; 

exhaust valve closure timing in-cylinder pressure 
calculation means for calculating an in-cylinder pressure at 
the exhaust valve closure timing; 
30 gas constant calculation means for calculating a gas 

constant corresponding to a change in a composition of 
exhaust gas, based on an air-fuel mixture ratio; 
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exhaust valve closure timing in-cylinder residual gas 
amount calculation means for calculating an in-cylinder 
residual gas amount at the exhaust valve closure timing, 
based on at least the in-cylinder temperature, the in- 
5 cylinder pressure, and the gas constant; 

valve overlap period blow-back gas amount calculation 
means for calculating a valve overlap period blow-back gas 
amount defined as a quantity of gas flow from one of intake 
and exhaust ports via a combustion chamber to the other port 
10 during a valve overlap period during which both of intake 
and exhaust valves are open together, the gas flow being 
created by a pressure difference between the intake and 
exhaust ports; and 

internal EGR amount calculation means for calculating an 
15 internal exhaust gas recirculation amount based on the 

exhaust valve closure timing in-cylinder residual gas amount 
and the valve overlap period blow-back gas amount. 

25. A method of estimating an internal exhaust gas 
20 recirculation amount of an internal combustion engine 

enabling the internal exhaust gas recirculation amount to be 

varied by changing a valve overlap, the method comprising: 
calculating an in-cylinder temperature at an exhaust 

valve closure timing based on an exhaust temperature; 
25 calculating an in-cylinder pressure at the exhaust valve 

closure timing based on an exhaust pressure; 

calculating a gas constant corresponding to a change in a 

composition of exhaust gas, based on an air-fuel mixture 

ratio ; 

30 calculating an in-cylinder residual gas amount at the 

exhaust valve closure timing, based on at least the in- 
cylinder temperature, the in-cylinder pressure, and the gas 
constant ; 
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detecting an engine speed; 

calculating a valve overlap period integrated effective 
area for an opening area of the intake valve and an opening 
area of the exhaust valve during the valve overlap period; 
5 determining the presence or absence of a supercharging 

phenomenon regarding the gas flow between the intake and 
exhaust ports via the combustion chamber during the valve 
overlap period, the gas flow being created by the pressure 
difference between the intake and exhaust ports, and the 

10 presence or absence of a choking phenomenon regarding the 
gas flow between the intake and exhaust ports via the 
combustion chamber during the valve overlap period; 

calculating a valve overlap period blow-back gas amount, 
which is defined as a quantity of gas flow from one of 

15 intake and exhaust ports via a combustion chamber to the 

other port during a valve overlap period during which both 
of intake and exhaust valves are open together, and is based 
on the engine speed, the valve overlap period integrated 
effective area, and the supercharging-and-choking phenomena 

20 decision result regarding the presence or absence of the 

supercharging phenomenon, and the presence or absence of the 
choking phenomenon; and 

calculating the internal exhaust gas recirculation amount 
by adding the valve overlap period blow- back gas amount to 

25 the exhaust valve closure timing in-cylinder residual gas 
amount . 

26. The method as claimed in claim 25, wherein: 

the valve overlap period integrated effective area is 
30 calculated based on the valve overlap. 



27. The method as claimed in claim 25, wherein: 
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the valve overlap period integrated effective area is 
calculated by integrating an effective opening area per unit 
time from a point of time corresponding to a beginning-of - 
valve-overlap crank angle to a point of time corresponding 
5 to an end-of -valve -overlap crank angle, the effective 

opening area being set to a smaller one of the opening area 
of the intake valve and the opening area of exhaust valve. 

28. The method as claimed in claim 25, wherein: 
the valve overlap period blow-back gas amount is 

calculated based on the engine speed, the valve overlap 
period integrated effective area, and a temporary valve 
overlap period blow-back gas amount determined based on the 
supercharging-and-choking phenomena decision result, from a 
predetermined expression: 

MRESOL = (MRESOLtm P xASUMOLx60)/(NRPMx360) 
where MRESOLtmp denotes the temporary valve overlap period 
blow-back gas amount, ASUMOL denotes the valve overlap 
period integrated effective area, and NRPM denotes the 
engine speed. 

29. The method as claimed in claim 28, wherein: 
an unsupercharged and unchoked state, temporary valve 

overlap period blow-back gas amount is calculated from a 
25 predetermined expression: 

MRESOLtmp = 4l xPEVCxMRSOLDxMRSOLP 
where PEVC denotes the in- cylinder pressure, MRSOLD denotes 
a density term and MRSOLP denotes a pressure difference term, 
and the density term MRSOLD and the pressure difference term 
30 MRSOLP are represented by predetermined expressions: 
MRSOLD = SQRT{l/(REXxTEVC) } 
MRSOLP = SQRT [ SHEATR/ ( SHEATR- 1 ) 

x{PINBYEX (2/SHEATR) -PINBYEX < (sheatr + i) /sheatr) } j 
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20 
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where REX denotes the gas constant, TEVC denotes the in- 
cylinder temperature, SQRT denotes a coefficient for the gas 
constant REX and the in- cylinder temperature TEVC, SHEATR 
denotes a ratio of specific heat of exhaust gas, 
5 corresponding to a change in a composition of exhaust gas, 
obtained when the supercharging phenomenon is absent, and 
the ratio of specific heat of exhaust gas is replaced with a 
ratio of specific heat of air-fuel mixture, obtained when 
the supercharging phenomenon is present, and PINBYEX denotes 
10 a pressure ratio of intake pressure to exhaust pressure. 

30. The method as claimed in claim 29, wherein: 

an unsupercharged and choked state, temporary valve 
overlap period blow-back gas amount is calculated from a 
15 predetermined expression: 

MRESOLtmp = PEVCxMRSOLDxMRSOLPC 
where PEVC denotes the in-cylinder pressure, MRSOLD denotes 
the density term and MRSOLPC denotes a pressure difference 
term, and the density term MRSOLD and the pressure 
20 difference term MRSOLPC are represented by predetermined 
expressions : 

MRSOLD = SQRT{l/(REXxTEVC)} 

MRSOLPC = SQRT[SHEATRx{2/(SHEATR+l) } {(SHEATR+1)/(SHEATR " 1)} ] 
where REX denotes the gas constant, TEVC denotes the in- 
25 cylinder temperature, SQRT denotes the coefficient for the 
gas constant REX and the in-cylinder temperature TEVC, and 
SHEATR denotes the ratio of specific heat. 

31. The method as claimed in claim 30, wherein: 

30 a supercharged and unchoked state, temporary valve 

overlap period blow-back gas amount is calculated from a 
predetermined expression: 
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MRESOLtmp = -V2xJ xPINxMRSOLPT 

\RINxTIN 

where PIN denotes the intake pressure, RIN denotes a gas 
constant of intake air, TIN denotes a Kelvin absolute 
temperature of the intake air, and MRSOLPT denotes a 
5 pressure difference term represented by a predetermined 
expression : 

MRSOLPT = SQRT [ SHEATR/ ( SHEATR- 1 ) 

x{PINBYEX ( - 2/SHEATR) -PINBYEX ( -<sheatr + i)/sheatr) } j 
where SQRT denotes the coefficient, SHEATR denotes the ratio 
10 of specific heat, and PINBYEX denotes the pressure ratio. 

32. The method as claimed in claim 30, wherein: 

a supercharged and choked state, temporary valve overlap 
period blow-back gas amount is calculated from a 
15 predetermined expression: 

MRESOLtmp = -J xPINxMRSOLPC 

V RINxTIN 

where PIN denotes the intake pressure, RIN denotes the gas 
constant of the intake air, TIN denotes the Kelvin absolute 
temperature of the intake air, and MRSOLPC denotes a 
20 pressure difference term represented by the predetermined 
expression : 

MRSOLPC = SQRT[SHEATRx{2/(SHEATR+l)} {(SHEATR+1)/(SHEATR " 1)} ] 
where SQRT denotes the coefficient and SHEATR denotes the 
ratio of specific heat. 

25 



